(N 

o 

(N 

u _ 

< 






cd 



I 

O 

o 



> 

m 
rn 

rn 

o 



X 



Topo-electronic transitions in Sb(lll) nanofilm: the interplay between quantum 
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When the dimension of a sohd structure is reduced, there will be two emerging effects, quantum 
confinement and surface effect, which dominate at nanoscale. Based on first-principles calcula- 
tions, we demonstrate that due to an intriguing interplay between these two dominating effects, 
the topological and electronic (topo-electronic) properties of Sb (111) nanofilms undergo a series of 
transitions as a function of the reducing film thickness: transforming from a topological semimetal 
to a topological insulator at 7.8 nm (22 bilayer), then to a quantum spin hall (QSH) phase at 2.7 nm 
(8 bilayer), and finally to a normal (topological trivial) semiconductor at 1.0 nm (3 bilayer). Our 
theoretical findings for the first time identify the existence of the QSH in the Sb (111) nanofilms 
within a narrow range of thickness and suggest that the Sb (111) nanofilms provide an ideal test 
bed for experimental study of topo-electronic phase transitions. 
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Thanks to the development of modern epitaxial-growth 
technique, the state-of-the-art electronic devices are often 
made of thin films with good control of film crystallinity 
and thickness with single atomic layer precision. When 
the film thickness approaches nanoscale, two emerging ef- 
fects become important. On one hand, the confinement 
in the z-direction normal to the film surface significantly 
alters the behavior of the itinerant electrons in the film by 
quantizing the electron wave, which manifests not only in 
film structural \V\ and mechanical stability [2], but also 
in electronic phases, such as magnetism j3|] and super- 
conductivity Q. On the other hand, the termination of 
a material with a surface leads to a change of the elec- 
tronic band structure differing from the bulk, with sur- 
face states formed only at the atomic layers closest to 
the surface [5|. Due to the large surface- volume ratio of 
a thin film, the surface effects may dominate the overall 
thin film electronic properties. 

Recently, a unique type of surface states is identified in 
materials with strong spin-orbit coupling (SOC), which is 
expected to serve as the dissipationless conducting chan- 
nels in spintronic devices [6|, |7|. Taking the form of a 
spin-resolved Dirac cone lying within the bulk band gap, 
the so-called helical surface states are theoretically at- 
tributed to the nontrivial Z2 topology of the bulk va- 
lence bands, which has in turn triggered a vast number 
of studies on these topologically nontrivial materials, i.e. 
the topological insulator (TI) in three dimensions (3D) 
and the QSH systems in 2D |i-[lll- 

An interesting surface coupling effect is observed in TI 
nanofilms 113, [l3. In such films, the helical surface states 



located on the top and bottom surfaces are sufficiently 
close to each other in space, so that the coupling be- 
tween the two surfaces becomes noticeable. The surface 
coupling effect opens a gap at the Dirac point on both 
surfaces and may even lead to topological transitions of 
the film [13, [H. 

One significant difference between the quantum con- 
finement and surface coupling is that they exhibit differ- 
ent scaling laws as a function of the film thickness. Fol- 
lowing the textbook description for quantum particles in 
a box, we have: 
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wliere AEb is the energy shift of a bulk state due to the 
quantum confinement, m* and L stand for the effective 
mass of charge carrier and the film thickness, respectively. 
The wavefunction of surface states takes the form of 
'^ss — u(r)e''^^^e**''''''e~^^, which modifies the general 
Bloch wavefunction in a periodic solid by replacing k^ 
with a complex number A. Tlie real and imaginary part of 
A are the decay constant and oscillating wave vector, re- 
spectively. The surface-coupling-induced splitting AEs 
is determined by {'i'ss \ '^'ss)^ where t and b stand for 
the top and bottom surfaces, respectively. Assuming the 
two surfaces are equivalent and by smoothing out u(r) as 
a constant, we have 
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When the thickness of a thin film is decreased, the 
quantum confinement with the "fiatter" power-law scal- 
ing is expected to show up first, modifying the bulk 
bands. The surface coupling will dominate in the 2D 
limit, because of the "steeper" exponential scaling. A 



joint action of the quantum confinement and surface cou- 
pling effect is likely to spawn rich physics in a nanofilm, 
especially when nontrivial topology is also involved. 

In this Letter, using first-principles calculations, we re- 
port our finding of a series of topo-electronic transitions 
in Sb (111) nanofilms triggered by the interplay of these 
two effects. As the film thickness reduces, the quantum 
confinement first opens a gap in the bulk band to induce 
a transition from a topological semimetal to a topologi- 
cal insulator. Then, the surface coupling effect kicks in, 
opening a gap in the surface bands giving rise to a QSH 
phase within a narrow range of thickness. The film finally 
degrades into an ordinary semiconductor in the 2D limit, 
because of the surface-coupling induced level crossings. 

Antimony is one of the main group elements with 
strong SOC. It crystallizes in a rhombohedral struc- 
ture as shown in Fig. la, which is typical among the 
group- V solids. Each atom has three equidistant nearest- 
neighbour (NN) atoms and three equidistant next-NNs 
slightly farther away. The atoms form a bilayered (BL) 
structure along the (111) direction (Fig. lb). Each bi- 
layer has a puckered honeycomb lattice with two atoms 
per unit cell, corresponding to the two NN atoms. A Sb 
(111) film consists of several BLs stacking up in an ABC 
sequence along the (111) direction. 
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thickness regime (1.3-2.7 nm), that is unknown before. 

Our calculations are performed with the density func- 
tional theory (DFT) using the plane wave basis, as imple- 
mented in the ABINIT package UM . We employ the local 
density approximation (LDA) (20[ and the Hartwigsen- 
Goedecker-Hutter pseudopotential [21|, which is gener- 
ated on the basis of a fully relativistic all-electron cal- 
culation and have been tested to be accurate for heavy 
elements like Sb. The spin-orbit coupling is incorporated 
in the self-consistent calculations as described in 2,2;]. 
The reliability of the standard DFT-I-LDA calculations 
on group-V semimetals including bulk Sb has been care- 
fully examined by comparing with experimental data for 
the density of states, number of free carriers and Fermi 
surface [23|. The numerical accuracy has been shown 
in strict control, with unexpected accuracy even for the 
Fermi surface, which is sensitive to the energy gap. 

Our calculated structural parameters of bulk Sb are 
shown in Fig. la and lb, in good agreement with previ- 
ous calculations [2j] . To model the thin film, a supercell 
of slab is used with periodic boundary conditions in all 
three dimensions and a 10 A vacuum layer between the 
slabs to eliminate the inter-slab interaction. The thin 
film has slightly relaxed structural parameters relative 
to the bulk values. All the atomic positions are fully re- 
laxed as the film thickness changes. A plane wave cut off 
of 28 Ry and a F-centered k-point mash of 8x8x1 are 
used. 
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FIG. 1: (Color online) (1) The hexagonal unit cell of single 
crystal Sb. (b) The top view of the Sb lattice. The colors and 
sizes of Sb atoms indicate different atomic layers. 




Bulk Sb has been predicted to be an intrinsic topo- 
logical semimetal [16|, which means that despite a neg- 
ative indirect band gap, the valence bands of bulk Sb 
are topologically nontrivial as those of a TI. The heli- 
cal surface bands in a 20-BL Sb film have been exper- 
imentally observed using the angle-resolved photoemis- 
sion spectroscopy (ARPES) [l3l- First-principles calcula- 
tions have also predicted a surface-coupling-induced gap 
in the 4-BL and 5-BL films [18|. However, a systematic 
study on all possible topo-electronic transitions in a Sb 
(111) nanofilm is still missing. Here, we map out its com- 
plete topo-electronic phases as a function of thickness, 
including the discovery of the QSH phase in a narrow 



FIG. 2: (Color online) The band structure of Sb (111) films 
with different thickness, (a) 26-BL, topological semimetal 
with a negative bulk gap; (b) 11-BL, TI phase; (c) 4-BL, QSH 
phase ; (d) 2-BL, ordinary semiconductor. The Fermi level is 
set to zero. For the 4-BL and 2-BL films, parity information 
at the r and M points is shown to distinguish their topological 
difference. 



In Fig. 2a, we plot the band structure of a 26-BL (9.4 
nm) film, which can be well understood as the bulk band 
projections (shaded regions) plus the surface bands (red 
curves) in the middle of the band gap. The conduction 
band edge and the valence band edge overlap at the Fermi 



level, showing the typical feature of a semimetal. A pair 
of surface bands cross at the F-point, forming a Dirac 
cone. The calculated band structure is in good agree- 
ment with recent ARPES data and earlier calculations 
[TA Il8| , which is consistent with that fact that bulk Sb 
is a topological semimetal. 

We now reveal how the quantum confinement and sur- 
face coupling effect modify the band structure when the 
film thickness decreases. For a semimetal, the dominat- 
ing effect of quantum confinement is to enlarge the en- 
ergy gap between the bulk conduction bands and valence 
bands due to the opposite signs of m* for electrons and 
holes. Two kinds of bulk gaps are traced in Fig. 3 as 
a function of the film thickness. One is the direct gap 
at the F-point; the other is the indirect gap between the 
conduction band minimum and the valence band maxi- 
mum. Both gaps are measured with respect to the bulk 
states excluding the two surface bands. 

The increase of both the direct (squares) and indirect 
gaps (diamonds) roughly follows the 1/L^ scaling as ex- 
pected, where L is the thickness. An important transi- 
tion occurs at 22-BL, where the indirect gap reverses its 
sign from negative to positive, signifying a semimetal- 
to-semiconductor transition, while the direct gap keeps 
open. Since the transition does not involve any level 
switching between the conduction and valence bands, the 
topology of the film does not change. Therefore, at 22- 
BL, the film transforms from a topological semimetal to 
a TI driven by the quantum confinement of bulk elec- 
tronic states. A typical band structure of the Sb TI film 
is plotted in Fig. 2b. 
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FIG. 3: (Color online) The direct (F point) and indirect bulk 
gap and surface splitting as a function of the film thickness. 
Inset :the enlarged plot of the surface band splitting from 3-BL 
to 8-BL. 

In Fig. 3, we have also plotted the splitting (dots) of 
the Dirac point formed by the two surface bands, which 
refiects the surface coupling effect. The Dirac point split- 
ting becomes noticeable only at a thickness much smaller 
than the semimetal-to-semiconductor transition point. 



This is because the exponential scaling of the surface 
coupling decays much faster than the power-law decay of 
quantum confinement. It is interesting to note that the 
splitting curve is not monotonic: first increases slowly, 
reaching the maximum at 4-BL; then drops to a minimum 
at 3-BL; finally increases rapidly from 3-BL to 1-BL (Fig. 
3) . A magnified view of the turn-around from 8-BL to 3- 
BL is shown as the inset of Fig. 3. Such an anomal y h as 
also been observed in earlier first-principles study [l^ , 
which attributed this behavior to the inversion-symmetry 
breaking. In our calculation, however, the inversion sym- 
metry is preserved for all the films. Therefore, some other 
mechanisms should be considered. 

Note that when deriving Eq. (2), we arbitrarily set 
the periodic Bloch function u(r) to be a constant. If the 
details in u(r) are retained, Eq. (2) should be rewritten 
as: 
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Here, we still assume that the two surfaces are equiva- 
lent and thus the in-plane components can be integrated 
out. Simply because of the rapid oscillating nature of the 
Bloch function, the additional integral in Eq. (3) has to 
be essentially non-monotonic as a function of L, varying 
not only in the magnitude but also in the phase. 

Besides a minor turnover as shown in Fig. 4, the mag- 
nitude modification does not override the overall expo- 
nential trend. The phase modification, however, may 
result in more significant changes in terms of the band 
topology [a,01- Specifically, for the Dirac point, u(r) can 
be chosen to be real, because it locates at the F-point. 
The phase, which is the sign in this specific case, of AEs 
in principle determines the relative positions of the bond- 
ing and anti-bonding levels formed by the coupled top 
and bottom surface states. It is known that a switch- 
ing between an occupied level and an unoccupied level 
with the opposite parity is going to change the overall 
topology of the occupied bands [la| . 

If we view surface-coupling-induced splitting at the 
Dirac point as the energy gap of a 2D system, the 
turnover of the surface splitting shown in Fig. 4 is analo- 
gous to a gap closing-reopening process, which is the pre- 
cursor of topological transitions [a, Q • To demonstrate 
this idea, we have calculated the Z2 topological invariant 
of Sb (111) film from 1-BL to 7-BL, which serves as the 
"order parameter" to differentiate the topologically triv- 
ial and nontrivial 2D phases. The calculation of the Z2 
invariant can be dramatically simplified by the so-called 
"parity method" [lg| , if the system is space inversion in- 
variant, as the case of the Sb (111) film. Accordingly, the 
Z2 number of the Sb film can be obtained from the wave- 
function parities at four time-reversal-invariant k-points, 
(Ki), one F and three M's, as: 
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FIG. 4: (Color online) The phase transition diagrams as a 
function of the film thickness from trivial semiconductor to 
topological semimetal. 



where ^ = ±1 is the parity eigenvalues and N is the num- 
ber of the occupied bands. The Z2 topological invariant 
v takes two values: v = 1 indicating a nontrivial phase 
and v — indicating a trivial phase. 



TABLE I: The total parity at the F and M points and the Z2 
number of Sb(lll) films with different thickness. 
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Table I presents the calculated Z2 topological invari- 
ants. Clearly, a phase transition occurs between 3-BL 
and 4-BL, which is in agreement with the speculation 
above. Between 1-BL and 3-BL, the films are topologi- 
cally trivial. Between 4-BL and 7-Bl, the films are topo- 
logically nontrivial, leading to a QSH phase. The calcu- 
lated surface band splitting of 4-BL and 5-BL is above 
SOmeV, providing a robust QSH gap at room tempera- 
ture. Typical band structures of the trivial semiconduc- 
tor phase and the QSH phase are plotted in Fig. 2c and 
2d, respectively. 



We should point out that the topology of the multi-BL 
Sb films cannot be interpreted as that of a stack of 1-BLs. 
Otherwise, all the Sb films would be topologically trivial 
because the 1-BL film is topologically trivial as found in 
Tab. L The outcome of the QSH phase is understandable 
after considering the level crossings induced by the inter- 
BL coupling as discussed for the Bi(lll) films 524|. 



In conclusion, by a systematic study of the band struc- 
ture and wavefunction parity as a function of thickness, 
we now have a complete understanding of the topo- 
electronic properties of Sb (111) films acrossing from the 
limit of 2D film to 3D bulk. As summarized in Fig. 4, in 
the "parameter space" of thickness, the topo-electronic 
phases can be divided into four distinct regimes. The ex- 
istence of such a rich topo-electronic phase diagram, un- 
derlied by a delicate interplay between the quantum con- 
finement and the surface effect, makes Sb (111) nanofilms 
an ideal test bed for experimental investigation of topo- 
electronic phase transitions. 
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